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We consider a pair of bosonic particles in a one-dimensional tight-binding periodic potential 
described by the Hubbard model with attractive or repulsive on-site interaction. We derive explicit 
analytic expressions for the two-particle states, which can be classified as (i) scattering states of 
asymptotically free particles, and (ii) interaction-bound dimer states. Our results provide a very 
transparent framework to understand the properties of interacting pairs of particles in a lattice. 
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Many-body quantum systems in spatially-periodic po- 
tentials were studied since the early days of quantum 
theory fl'] . Several simple but powerful models have been 
invented and successfully employed for the description of 
various condensed matter systems, including the Heisen- 
berg spin and Hubbard models. Lately, the Hubbard 
model has been attracting renewed attention prompted 
by the remarkable experimental progress in cooling and 
trapping bosonic and fermionic atoms in optical lattices 
0. The relevant parameters of this system can be con- 
trolled with high precision and can be tuned to imple- 
ment the Hubbard model with unprecedented accuracy 
i- 

Recently, Winkler et al. [4] have observed a remark- 
able Hubbard model phenomenon: A repulsive on-site 
interaction between a pair of bosonic atoms in an optical 
lattice can bind them together into an effective, dynam- 
ically stable "dimer" 0, [a, 0] ■ Here we revisit the prob- 
lem of two bosonic particles in a ID (one-dimensional) 
periodic potential described by the Hubbard model. We 
consider both cases of attractive and repulsive on-site in- 
teractions between the particles and derive simple, exact 
analytic solutions for the two-particle problem. One type 
of these solutions corresponds to the scattering states of 
asymptotically free particles, while the other describes 
the interaction-bound dimer states with energies below 
or above the energy band of the scattering states [1| . 

We note that the solution of the Schrodinger equa- 
tion for the two-particle problem in a lattice can be de- 
rived using the Green function approach Q, which 
although being more general, requires laborious calcula- 
tions often accompanied by approximations for obtaining 
explicit analytic results. In contrast, in what follows we 
give an extremely simple and straightforward derivation 
of the energies and exact analytic wavefunctions for the 
two-particles states in a lattice. 

The Hamiltonian of the problem is given by 



H 



-J 



bj) + 



U 



(1) 



where foj (bj) is the bosonic creation (annihilation) op- 
erator and fij = b^jbj the number operator for site j, 
J(> 0) is the tunnel couphng between adjacent sites, 
and U is the on-site interaction. A standard basis 



{ \nj)} for Hamiltonian ([T]) is composed of the eigenstates 
\nj) = --i=(&j)" I vac) of operator fij whose eigenvalues 
n = 0, 1, 2, . . . denote the number of bosonic particles at 
site j, and |vac) = |{0j}) is the vacuum state. 

For a single particle in the periodic potential, the last 
term on the right-hand side of Eq. ([T]) does not play a 
role. The Hubbard Hamiltonian then takes the form 



(2) 



where \xj) = \lj) denotes the state with a single parti- 
cle at position Xj = dj of the jth lattice site, with d the 
lattice constant. Expanding the single-particle state vec- 
tor as \^p) = J^ji^i^j) the stationary Schrodinger 
equation H'-^^ lip) — E^^^ \ip) reduces to the difference 
equation 

- J[i^{xj^i) + H^j+i)] = E^''> i:(x,) (3) 

for the wavefunction ip(xj). Taking ipi^) — ''Pqi^) = 
exp{iqx) immediately yields the well-known result 

- 2 Jcos(gd) ^,{xj) = M^3)^ (4) 

which implies that the discrete plane waves ipqi^j) = 
exp{iqdj) are the eigenf unci ions of Hamiltonian ([2]), with 

the corresponding eigenenergies E^'^^ = -2Jcosiqd), 
forming a Bloch band of width 4 J [l| . The effective mass 
TO* of the particle with quasi-momentum q close to is 
then given by the usual expression 



d^E, 



(1) 



dq'^ 



q=0 



2Jd^' 



(5) 



For two particles in a one-dimensional lattice, the Hub- 
bard Hamiltonian can be recast in the explicit form 



-J2^i\Xj){x3 + l\ + \Xj+i){x^\) 
j 

+ \yj+i){yj\) 



(6) 
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while the two-particle state vector can be expanded in 
terms of the non-symmetrized basis { \xj, yj>)} as 



(7) 



where yj = dj. Clearly, the standard (symmetrized) 
bosonic basis is related to the non-symmetrized basis 



via the transformation \2j) — \xj,yj) and 



7f(kj,yi') + \yj,Xj')) (j ^ /). The stationary 

Schrodinger equation H'^^^ j'l') = E'-^^ I'i') is now equiva- 
lent to the recurrence relation 

-J[i'ixj^i,y,,) + -9{x,+i,y,,) 
+^ix,,y,,^i) + -9{x,,y,,+i)] 

+C/%,*(x,,y,0 = i?(2)vI/(x„y,0. (8) 

Define the center of mass R — ^{x + y) and relative 
r = X — y coordinates. In terms of the new variables, the 
two-particle wavefunction can be separated as 



(9) 



where the relative coordinate wavefunction ipK{r) de- 
pends now on the center-of-mass quasi-momentum K e 
[— Tr/d, Tr/d] as a continuous parameter. Equation ([8]) 
then yields 

- JK[ipKiri-i) +ipKiri+i)] + USro-ipKiri) 

= e'^^ ^JjKin), (10) 

with Jk = 2J cos{Kd/2) and — di {i — j — j')- 

(i) Consider first the scattering solutions for a pair of 
asymptotically free particles [1, Q. Since the action of 
the short-range scattering potential U5ro amounts to a 
unitary phase-shift (see below) , the spectrum of such so- 
lutions is given by the sum of the spectra for two free 
particles x and y with momenta Qx = K/2 + k and 
Qy = K/2 - k: 



= -4Jcos(i\:d/2)cos(fcd) 

= -2JKC0s{kd). (11) 



Obviously, this result also follows directly from Eq. (|T0|) . 
with U = 0, upon substitution of the plane waves 
i'Kir) = "ipK ki^) ~ exp(±iA:r) For a given value 
of the center-of-mass momentum K, and thereby Jk, 

(2) (2) 

the lowest Ej^Q = —2Jk and highest E)^'^ = 2Jk en- 
ergy states of a pair of asymptotically free particles cor- 
respond, respectively, to the relative momenta fc 
and k ^ n/d. The energy spectrum of Eq. (fTT|) . and the 
corresponding density of states 



L dk L 



1 



271 dE 27rd ^[Ucos{Kd/2)Y - E^ 
with L a quantization length, are shown in Fig. [TJ 



, (12) 



We can now derive an explicit expression for the wave- 

(2) 

function of the scattering states. Substituting E)^\. of 
PTjl into Eq. pUj) . setting ipK.kiO) — C and using the 
bosonic symmetry ipK,kir) — '4'K,k{—'r), after little alge- 
bra we obtain 

ipK.k{ri) = Ccos{kri) + C^^^TT^— ^ sin(/c|rj|). 

In terms of the plane waves undergoing the scattering 
phase shift 5K,k, the wavefunction reads 



i^K,k{ri) 
ta.n{5K.k) 



U csc{kd) 



(13) 



2J 



K 



In the limit of [/ — > 0, this reduces to the trivial solution 
''pK.ki'fi) = cos{kri) for two non-interacting bosons, while 
in the opposite limit of U/Jk ±oo, Eq. yields 
the fermionized solution i^K^kifi) = sin(fc|ri|). At the 
edges of the scattering band fc ^ 0, Tr/d, we can define a 
generalized ID scattering length ax via Q 



aji = — lim 



86 



K.k 



2dJ 



K 



fc^o dk 



U 



(14) 



which is thus positive for [/ < and negative for J7 > 0, 
as can be intuitively understood on the basis of pertur- 
bation theory Q. Thus, in the case of L/ < 0, a pair of 
co-localised particles (forming a dimer discussed below) 
has an energy below that of a pair of unbound particles, 
causing therefore an effective repulsion between the un- 
bound particles by pushing their energy upwards. Con- 
versely, the energy of the repulsively bound pair [/ > is 
above the energy of the unbound pair of particles, which 
lowers their energy causing an effective attraction. 

(ii) The on-site interaction [/ 7^ can bind the two 
bosonic particles together into a close dimer 0, d, 0| 
whose energy is above [U > 0) or below {U < 0) 
the continuum (fTTjl of scattering states. We now de- 
rive simple analytic solutions for the bound dimer states. 
For \K\ = n/d {Jk = 0), Eq. ^ immediately yields 



E^^l = U and 



i^./din ^ 0) = 0. 



(15) 



For K G (— Tr/d, 7r/d), using the exponential ansatz 
where C is a normalization constant. 



from Eq. (flOl) we have 



-2JKaK + U = eP, 

, \i\ + l 



a 



K 

p(2) 



K ' 



(16a) 
(16b) 



where the unknowns arc E'^' and UK- The solutions for 
aK are given by 



aK=UK± JUf^ + 1, Uk 



U 
2J^' 
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(2) 

FIG. 1: Energies E)^ versus the center-of-mass momentum K 
for a pair of bosons in a ID lattice described by the Hubbard 
model. The continuum spectrum corresponds to energies (|lip 
of the scattering states, with the shading proportional to the 
density of states (fT2|) . The (blue) line below and the (red) line 
above the scattering band are, respectively, the energies of the 
attractively-bound dimer with U = —5J and the repulsively- 
bound dimer with U = 5J. Their relative coordinate wave- 
functions (|17b|) and (|18bl) at \K\ = 0, 7r/2d, n/d are shown on 
the bottom and the top panels. 



A normalizable bound state corresponds to \aK\ < 1, 
which exists for any U ^ 0. Thus, for attractive interac- 
tion, U < 0, we have 



aK 



1-\Uk\>Q. 



with which for the energy and normaUzed relative coor- 
dinate wavefunction we obtain 



E 



(2) 
K 



(17a) 



U'ji + l-\UK\) . (17b) 



The lowest energy of an attractively-bound dimer is 
-VC/^ + 16J2 attained aX K = while the high- 



E 



(2) 



est energy is E)^l^ = —\U\ = U aX, K = ±-K/d. We can 
define the binding energy for a dimer with momentum K 
as the difference between the energy of the bound pair 
E)^' and the energy at the bottom of the scattering band 



for an unbound pair E^^\^ ^ -2Jk (see Fig. [H central 
panel) , 



C/2 + 4 J|. < 0. 



For repulsive interaction, J7 > 0, we have 



OlK = Uk 



1 < 0, 



and the dimer energy and wavefunction are given by 



E 



(2) 
K 




(18a) 
(18b) 



The lowest and highest energies of a repulsively-bound 



dimer are given by E^j'^ 

VU'^ + 16J^, attained, respectively, at K ~ ±.TT/d and 
K ~ Q. The dimer binding energy, counted from the top 
of the scattering band E\.'^ = 2Jx (see Fig. [T]), is now 
positive, given by 



\U\ 



U and £'^- 



(2) 



E 



(2) 
K 



^(2) 



2Jk > 0. 



Clearly, for a given value of momentum K of the dimer, 
the stronger the on-site interaction \U\, the smaller the 
extent of the dimer wavefunction. In the case of re- 
pulsive interaction U > {uk < 0), the sign of wave- 
function (jl8b[) alternates between the neighboring sites 
i. Remarkably, when \K\ — Tr/d, and thereby Jk — 
0, the dimer wavefunction is completely localized at 
r,; — for any U ^ 0. To illustrate the foregoing 
discussion, in Fig. [1] we show the dispersion relations 
(fT7a|l and and the wavefunctions (jlTbp and (|18b|) 

at \K\ = 0, 7r/2(i, Tr/d for attractively- and repulsively- 
bound dimers. 

The dimer effective mass M* for small K is given by 



M* = 



d^E 



J K=0 



where the upper sign corresponds to attractive interac- 
tion U < 0, leading to positive M*; while the lower 
sign stands for repulsive interaction [/ > 0, for which 
the effective mass M* is negative. For weak interaction 
\U\ < J, we have M* ~ ±h^/{(Pj) = ±2m*, i.e., twice 
the single particle effective mass m* of Eq. ([5]). On the 
other hand, for strong interaction \U\ 3> J, we obtain 
M* ~ ;iV(2dV(2)), where J^^^ = -2P/U is the effec- 
tive tunnelling rate of the dimer between the neighboring 
lattice sites ^, iS, lQ, lL] • Now the dimer effective mass is 
large due to its slow tunnehng ^ J, with the sign 

of J^^-' determining the sign of M* . In this limit the 
energies (II Tap and (|18ap can be approximated as 



E^^^ ~{U- 2J(2)) - 2J(2) cos(Kd), 



(20) 
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where the first term on the right-hand side represents 
the dimer "internal energy" , while the second term is the 
kinetic energy of a dimer with quasimomentum K. 

To summarize, we have considered a pair of bosonic 
particles in a ID tight-binding periodic potential de- 
scribed by the Bose-Hubbard model and derived ex- 
act analytic solutions for the two-particle problem. 



These solutions describe the scattering states as well as 
interaction-bound dimer states of the two particles. We 
have obtained simple explicit expressions for the energies 
and wavefunctions of the two-particle states and outlined 
their properties. Our results are relevant to the recent ex- 
perimental 4] and theoretical 0, Q studies of strongly 
interacting pairs of atoms in an optical lattice. 
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